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ABSTRACT 

In the absence of constraints from the binary companion or supernova remnant, the 
standard method for estimating pulsar ages is to infer an age from the rate of spin-down. 
While the generic spin-down age may give realistic estimates for normal pulsars, it can 
fail for pulsars with very short periods. Details of the spin-up process during the low 
mass X-ray binary (LMXB) phase pose additional constraints on the period (P) and 
spin-down rates (P) that may consequently affect the age estimate. Here, we propose a 
new recipe to estimate millisecond pulsar (MSP) ages that parametrically incorporates 
constraints arising from binary evolution and limiting physics. We show that the stan- 
dard method can be improved by this approach to achieve age estimates closer to the 
true age while the standard spin-down age may overestimate or underestimate the age 
of the pulsar by more than a factor of ~10 in the millisecond regime. 

We use this approach to analyze the population on a broader scale. For instance, 
in order to understand the dominant energy loss mechanism after the onset of radio 
emission, we test for a range of plausible braking indices. We find that a braking index 
of n=3 is consistent with the observed MSP population. We demonstrate the existence 
and quantify the potential contributions of two main sources of age corruption: the pre- 
viously known "age bias" due to secular acceleration and "age contamination" driven 
by sub-Eddington progenitor accretion rates. We explicitly show that descendants of 
LMXBs that have accreted at very low rates (?h <C M^dd) will exhibit ages that appear 
older than the age of the Galaxy. We further elaborate on this technique, the impli- 
cations and potential solutions it offers regarding MSP evolution, the underlying age 
distribution and the post-accretion energy loss mechanism. 

Subject headings: pulsars: general — stars: neutron — stars: statistics — X-rays: 
binaries 



^FuU resolution color figures and movies available at URL: http://www.kiziltan.org/research/MSP/ages.html 
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Introduction 



An accurate determination of pulsar ages plays a critical role in our understanding of advanced 
stages of stellar evolution, supernova explosions and remnants, white dwarf (WD) atmospheres and 
cooling models, binary evolution, planet formation around compact objects, and pulsar evolution 
in general. 

Typically, pulsar ages are estimated by calculating the amount of energy lost during their 
spin-down. Consequently, the spin-down age of a pulsar can be formulated as: 



(n-l) P 



1- ^ 



P 



(1) 



where the period (P) and the spin-down rate (P) are the two main observables acquired by pulsar 
timing measurements. In the standard approach, the unknown initial spin period (Pq) of the pulsar 
is assumed to be much smaller (Pq <C P) than the observed period. The dominant energy loss 
mechanism is analytically captured by the braking index, which is n = 3 for pure dipole radiation 
and is implicitly adapted by the characteristic age r^. The age of a pulsar can then be conveniently 
approximated by its characteristic age, where: 



P 
2~P 



for Po<^P, n = 3. 



(2) 



While this approach may give reliable estimates for some normal pulsars (e.g. Crab Pulsar- 
PSR B0531+21: Tc ^ 12 40 yr whereas the age from the supernova remnant (SNR) tsnr ~ 955 yr; 



Wvckoff &: Murravlll977l ). it should be kept in mind that characteristic ages for some other pulsars 
will suffer, dramat ically in some cases (e.g. PSR J0205+6449: Tc ~ 5370 yr, tsnb. ~ 820 yr; 
Murrav et al.ll2002l ). because of the assumptions that render the standard approach less accurate, 
especially in the millisecond regime. Therefore, it is useful to develop a more comprehensive and 
detailed framework by which we can quantitatively understand MSP spin evolution. 

In our previous work, we have set up a base by which we demonstrated that the proper 
i nclusion of evolutioiiary co nstraints alone gives a deeper insight into the subsequent spin evolution 
(|Kiziltan fc Thorsettl[2009l ). 



In this paper, we propose a recipe to estimate pulsar ages that parametrically incorporates 
additional evolutionary and physical constraints. We show that the combined effect of a possible 
spin-up process during the LMXB phase and a maximum spin period due to the limiting centrifugal 
forces imparts meaningful constraints on the joint period (P) and spin-down (P) values that MSPs 
can attain, which ultimately are used to estimate their ages. We detail the contribution this new 
approach offers to our understanding of MSP evolution and elaborate on the ramifications on several 
related problems such as the dominant energy loss mechanism and braking indices of millisecond 
radio pulsars, WD atmospheres and cooling models, the underlying age distribution, the enigma of 
MSPs that appear older than the galaxy they reside in, and the sources of MSP age corruption. 
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The Spin Up Process 



The period and spin-down (P — P) relation of a particular MSP at the epoch when it turns 
on its radio emission can be scaled as 



P oc rhPr 



4/3 







(3) 



as the neutron star can ultimately be s pun-up no further tha n the spin period delineated by the 
Keplerian velocity at the Alfven radius ([Ghosh &: Lamblll992l ). We will use the least conservative 
upper bound of this scaling factor as an upper limit for the region where MSPs may be re-born. 
A single birth line instead would imply that these MSPs have accreted at near-Eddington rates 
(m ~ MEdd), which appears unreasonable, at least for a significant fraction of the observed MSP 
population as indicated by the paucity of sources near the spin- up line (see i i4.1l for discussion). 

We also now know that it is more than l ikely that the observed ma jority of MSPs do not initially 
re-appear in the vicinity of the spin-up line (jKiziltan &: Thorsettll2009l ). The region where MSPs are 
re-born as radio sources on the P — P diagram is strongly correlated with the dominant accretion 
rate (m) experienced during the last phases of LMXB evolution, which is poorly constrained. For 
the scope of this paper, using the spin-up line as a marginal upper boundary for the region where 
MSPs are re-born as radio sources (i.e. where the spin-down trajectories on the P — P plane 
start) will adequately account for possible nonlinear offsets due to other inherent uncertainties and 
assumptions made in Equation ( [3]) regarding the accretion geometry (streams of hot plasma flowing 
onto the neutron stars' polar caps instead of uniform spherical or wind accretion) or the opacity of 
the accreted material (sensitive to whether the companion has an H or He rich envelope attached 
to a CO or ONeMg core, see §: [3X2]) . 

On the other hand, the cumulative uncertainty of a presumed spin- up line t hat would offset o r 



1982 



Bhattacharva fc van den Heuvel 



1991 



tilt a particular re-birth line cannot be arbitrarily big (lArzoumanian et al.lll999l : 
Also, within the context of the standard recy cling scenario (jAlpar et al. I ll982l : 



Frank et al 



200211. 



Radhakrishnan &: Srinivasan 



, see for review), the "spin-up line" is merely an upper 
boundary below which MSPs are expected to be born rather than the line of culmination. To this 
day, there are no recycled pulsars observed above the spin-up line, except few in globular clusters 
(GCs) which have very uncertain evolutionary histories. Therefore, we will limit ourselves to Galac- 
tic MSPs whose spin-down history, orbital dynamics and Galactic kinematics remain unperturbed 
by gravitational encounters. 



3. Millisecond Pulsar Ages 

The standard approach to estimate pulsar ages in the absence of additional constraints from 
either a possible association to an SNR or a stellar companion, has been to use the characteristic 
age as a proxy to the true age. The main goal of our work is to better understand the non- 
trivial relationship between the physically important true age and the observationally accessible 
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characteristic age. We will refer to the time that has passed since the cessation of accretion as the 
"(true) age: rt" of a recycled pulsar. 



3.1. Alternative Methods to Estimate Pulsar Ages 

3.1.1. Ages From Kinematics and Supernova Remnants 

For some young pulsars that have reliable proper motion and distance measurements, a kine- 
matic age estimate can be made by tracing the pulsar's trajectory in the galactic gravitational 
potential. However, without a firmly established birth site, kinematic ages are at best an indirect 
means to constrain pulsar ages. 

For MSPs, that have no associated SNRs and have ages that are much longer than the orbital 
timescales in the Galactic potential, all kinematic age information has been lost. 



3.1.2. White Dwarf Cooling Ages 



After the discovery of optical emission from pulsar companions (jKulkarnil 119861 ) . WDs were 
soon realized as an alternative means to estimate the age of an MSP. Once active accretion has 
ceased in LMXBs with a recycled pulsar primary, the WD begins to subsequently radiate its internal 
heat after it burns off the remaining envelope. So, the beginning of WD cooling also marks the 
epoch when the spin-down starts for its companion. Therefore, in pr inciple, the WD cooling a ges 
are expected to be consistent with the ages of their MSP companions (jHansen &: Phinnevlll998al lbl). 



The basics of WD cooling models are potentially accessible to theoretical understanding be- 
cause of the simple thermal structure of the WD. The whole system is kept isothermal due to the 
efficient heat conduction of degenerate electrons. However, in practice, cooling ages remain difficult 
to estimate once realistic (and complex) effects of surface physics and stellar structure are included , 
leading to discrepancies and controversies in the interpretation of spec ifi c observations (IWoodlll992l : 
Sarna et al.lll998l . boool : ISchonberner et ahlboool : lAlthaus et aPboOlbl lal: Ivan Kerkwiik et ahlboosh . 



Possibly more promising than the original goal of using WD cooling ages to constrain the 
properties of MSPs, one might instead hope to use better constrained MSP ages (r, see §: 13. 3|) to 
understand WD atmospheres and cooling models. 
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3.2. Characteristic Ages: Idealized Pulsar Spin-down 



For millisecond radio pulsars, unbiasecj^ characteristics ages (r^) become upper limits to the 
true ages (r^) as the spin-down trajectories are truncated below the spin-up line. MSPs can be 
re-born on the spin-up line only if they accrete at an Eddington rate (m = MEdd) during recycling 
(see E quation ([3l)). Because of li kely sub-Eddington accretion rates experienced during the LMXB 



phase (IKiziltan &: Thorsettll2009l ). the majority of the spin-down trajectories start well below the 
spin-up line. A considerable fraction of millisecond radio pulsars (~30%) would even be expected 
to be born below the Hubble line (see §: 14.11 for discussion). In this approximation, Tc is derived by 
implicitly assuming pure dipole spin-down in the absence of other forms of additional torques and 



( Bildsten 


1998: 


Krolik 


1991 



are 

also assumed to be absent when Tc is inferred from the observed P — P values. Possible non- 
monotonic field decay before magnetic stability sets in may also contribute to the corruption of 
characteristic ages. 

O ne bias for which we can properly correct is the effect of secular acceleration, i.e., Shklovskii 
effect (jShklovskiil 119701 ). The observed P values include an additional apparent spin-down factor 
introduced because of the increasing projected distanc e between the puls ar and the solar system 
barycenter. This leads to a quadratic centrifugal term (jCamilo et al.lll994l ). 



A = 1.08 X 10" 



-18 
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m 



P ^ P" ^ 



p p 



(4) 



(5) 



where P and P' are the measured and unbiased spin-down rates for a pulsar at a distance D (in 
units of kpc) with a transverse velocity vt (in units of 100 km s~^). 

Figure [1] shows the observed MSP population and the extent of the bias introduced by secu- 
lar acceleration. MSPs that have a combination of relatively high transverse velocities and small 
distance measurements, have higher corresponding correction factors, e.g. for PSR J0034— 0534 
(D=0.98 kpc, Vt ~146.3 km s-^) and PSR B1257-M2 (D=0.77 kpc, vt ~350.6 km s'^) the corre- 
sponding correction factors are P jP' ~9.3 and 16.9 respectively. 



^For consistency, we designate unbiased values that are corrected for secular acceleration (i.e., Shklovskii effect) 
by adding " ' " to the parameter in lieu of referring them as "intrinsic" values. While the unbiased spin-down rates 
will represent the intrinsic values (i.e., P' — Pi), the unbiased characteristic age = P/2P' is neither the intrinsic 
nor the true age (i.e., 7^ n = rt, see 14. II for discussion). 
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3.3. A Realistic Age for Millisecond Pulsars (r) 

One of the reasons why characteristic age estimates become less reliable for MSPs, even after 
the observed spin-down rates are unbiased for secular acceleration, is due to the assumption that 
the birth periods are much smaller than their currently observed spin periods (Pq ^P), which fails 
for a considerable fraction of the population (see § I4.ip . H 



In fact, the predicted MSP spin-down age is proportional to the integrated (and normalized) 
spin-down path from Pq to P. The difference in the integrated trajectories will depend on the initial 
spin period Pq, which remains an unknown parameter in most cases. We know that the spin-down 
trajectory in the case when recycling and a constraining maximum spin limit are taken into account 
is more compressed than what the standard approach predicts. Therefore, with tighter upper limits 
on the true age, r will give an age estimate that is closer to r^. 

Given enough time, MSPs may reach a limiting equilibrium phase where they begin to lose 
angular momentum gained from accretion by shedding mass. A maximum spin limit beyond 
which a neutron star cannot be spun-up due to this continuous loss of excess angular momen- 



tum will truncate the spin-down trajecto ries vertically. Several authors (jHaensel &: Zdunik 



Friedman &: Ipser 



1992 



198S 



Cook et al.lll994l ) have constructed equilibrium sequences for neutron stars 



where a range of mass shedding periods are calculated for different equations of states including the 
effects of rapid rotation and large deviations from spherical symmetry. While the theoretical best- 
fit values range between Psh =1.28-1.32 ms for neutron stars with realistic confi gurations, for some 



extrem e cases they find that a limiting period of Psh ^0.85 ms may be plausible. IChakrabarty et al 
([2003) find evidence for a statistically significant upper limit at Pgh —1.32 ms. We therefore use a 
putative value of Psh ~1 and include it parametrically to perform calculations. 

The critical magnetic field Be will be the locus of points extending from the intersection 
of the diagonal spin-up line and the vertical mass shedding limit. MSPs with magnetic fields 
below Be can only be spun-up to the limiting mass shedding period. The critical magnetic field 
is Bc= (3.36^gg) X 10® G for P=l ms where the spin -up line is prescribed as -P = aP^^^ for 
a = (1.1 ± 0.5) X 10~^^s~^/3 (jArzoumanian et al.l[l999l ). 



We can formulate an age estimate that implements the constraints arising from recycling and 
mass shedding: 



t(B <Br 



P 



(n-l) P 
P 

(n-l) P 
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n-l 
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'Psh 



n-l 



(6) 
(7) 



^See a time-lapse movie for the true age evolution of millisecond pulsars at URL: 
http://www.kiziltan.org/researcli/MSP/ages.html 
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where is the mass shedding Hmit. We parametrically adopt the re-normahzed coefficient a = 
2.6tol X lO^s^/'^ and use it as a fiducial value. Although a inherently has numerous sources of 
uncertainty, the corresponding minimum post-accretion period that defi nes the spin-up line depend s 



on the total accreted mas s, but is insensitive to other parameters (see lPhinnev &: Kulkarnilll994l ) 



Arzoumanian et al.l ( 19991 ) find an upper limit of a = 1.6 x 10 ^^s to be empirically reasonable 
which we adopt to marginalize 3. 

At low magnetic fields (B< Be), the Alfven radii (r^ cc Bs ) shrink essentially down to the 
surface of the neutron star and angular momentum is transferred more efficiently and stably on 
longer timescales. Therefore the spin-down trajectories for neutron stars with lower magnetic fields 
are more likely to start closer to the vertical re-birth line (~Psft). 

The level of truncation of the spin-down trajectories can be more dramatic for neutron stars 
with very low magnetic fields and short periods (P< 3 ms). These MSPs with B< 10^ G are 
expected then to be born with periods very close to their observed ones (P~Po) and therefore are 
more likely to be much younger than they appear. 

In Figure [T] the sets of blue and red lines are MSP age (r) lines for braking indices n=3 
(red: dash) and 5 (blue: dash-dot). It is noteworthy to point out that for MSPs whose angular 
momentum and energy loss is dominated by multipole or quenched by gravitational wave radiation, 
the scaling factor in Equation ( [3]) can be different. The blue line in Figure [1] demonstrates the 
potential level of additional bias for the case in which gravitational wave radiation is the dominant 
mechanism for energy loss as opposed to braking due to pure dipole radiation. MSPs that lose 
energy via more efficient processes will traverse the spin-down path much faster, and consequently 
mimic older ages. As a result, the contribution of more efficient processes will exacerbate the age 
overestimate further. 

Figure [2] shows t he true age trend for MSPs . The synthetic population is produced with the 



method described in iKiziltan &: ThorsettI (|2009l ). where the MSP evolution is parameterized by 



the evolution functional £(D,M,R), in which D is the initial period (Pq) distribution of the 
progenitor population, M is the distribution of predominant accretion rates experienced during the 
latest phases before the onset of radio emission, and R is the galactic birth rate. The parameter 
space {D, M, R) is sampled by producing purely random synthetic populations. Then, only the 
input parameters that produce synthetic samples consistent with the observed MSP population are 
used for the construction of the underlying P — P demographics. For the filtration process, we 
use a relaxed mult i -dimensional Kolmogorov-Smi rnov (K-S) test as the cons istency criteria 



for the multi-layered Monte-Carlo integration scheme (jFasano &: Franceschinilll987l . and references 
therein) . 

We choose sub-samples that are uniformly selected from the underlying MSP population to 
produce Figures [2] and [3l While selection biases play a role in which MSPs are preferentially 
observed, the reflected age trend will remain unaffected. Therefore, the age trend of the underlying 
population both in Figures [2] and [3] is expected to be a realistic reflection of the true age. 
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4. Discussion 

4.1. Sources of Age Corruption: Bias and Contamination 

Older MSPs tend to appear younger if the spin-down rates (P) are not properly corrected for 
the contribution due to secular acceleration. This causes an upward "age bias" on the P — P plane. 
The Shklovskii effect is more pronounced for MSPs that are closer to the solar system barycenter 
and have higher transverse velocities. 

We can unbias our measurements and correct for this corruption by calculating the centrifugal 
term (Equations and ([5])) once both the distance and transverse velocity terms are accurately 
known. The proper motion measurem ents for MSPs with poorly constrained distances pose a stiffer 



challenge than previously predicted (jPeller et al.ll2009l ). To avoid underestimating the potential 



age bias, we introduced conservative random deviati ons to the observed distance and transverse 



velocity measurements by up to ±OAD and ±0.05^4 (jBrisken et al.ll200ll2003l : IPeller et al.ll2009l ) 
to produce the sample population in Figure O 

Figure [3)^a) reflects the potential age bias we estimate for the underlying population. In some 
cases, relatively old MSPs may even appear above the spin-up line if the observed spin-down rates 
are not properly unbiased. The age bias caused by the Shklovskii efl^ect will tend to push older 
MSPs upward by amounts proportional to their correction term [P/P') and hence make them 
appear younger. 

Throughout this work, we have excluded MSPs in GCs due to their complicated spin evolu- 
tion. The cluster's compact gravitational well, and to a lesser extent the larger cross section for 
gravitational interaction in GCs relative to the Galactic disk, efl'ectively perturb the spin evolution 
of these MSPs by injecting a cumulative corruption to their spin-down rates. Conversely, the spin- 
down trajectories of MSPs in the Galaxy will include less corrupted information about the initial 
formation. Even though the level of age bias for Galactic MSPs appears not as dramatic as the 
ones that are gravitationally perturbed in GCs, it still mixes the apparent age distribution quite 
efficiently (Figure Et^a) ) . 

Unbiasing will push old MSPs that appear younger to their corresponding age (r) lines. This 
correction process will recover the final positions of the true spin-down trajectories but will eff'ec- 
tually exacerbate the downward age contamination which was artificially diluted by secular accel- 
eration. The main source of the downward age contamination though will be the sub-Eddington 
progenitor accretion rates experienced during the LMXB phase. 

These MSPs are expected to be born with smaller spin-down rates and consequently traverse 
shorter trajectories. We predict that ~30% of MSP ages overestimate the true age by more than 
a factor of 2. Therefore, we argue that MSPs, which were presumed as the oldest sub-population, 
have a fiatter true age distribution than previously thought. 

While we can reverse the upward "age bias" , the downward "age contamination" on the other 
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hand, reflects upon the intrinsic spin-down rates and is real (Figures [3]^a) and (b)). Unless we 
have a unique leverage to accurately determine the period at birth (Pq) along with the respective 
progenitor accretion rate, the inferred MSP ages will remain to be strict overestimates (Figure [3)^b)). 

Figure 131 quantifies the level of bias and contamination of MSP ages before and after correcting 
for the contribution due to secular acceleration. Both the upward bias and downward contamination 
will simultaneously remain before the spin-down rates are properly corrected for the Shklovskii 
effect. We predict that 20% of the measured P — P values for MSPs (dotted line in Figured]) will 
overestimate the age by more than a factor of two, whereas about 10% will underestimate the age 
at the same level. Old and young MSPs are practically indistinguishable before the spin-down rates 
are properly unbiased. Some very young MSPs may appear even below the Hubble line whereas 
much older MSPs may be observed above the spin- up line (Figure [3]^a)). The asymmetric wings of 
the dotted line that extend in both directions in Figure [4] quantify this bidirectional corruption. 

The ages obtained from intrinsic P — P values will represent strict upper limits to the true 
age (tj = Tt <t' < T^). The truncated solid line in Figure H] shows at what level we overestimate 
the ages for the whole MSP population. After properly unbiasing the observed spin-down rates, 
we expect to overestimate the age of 30% of MSPs by more than a factor of two. Table 1 shows Tc 
and T for observed millisecond radio pulsars before and after unbiasing. Table 2 marginalizes the 
potential bias for an assumed Vt -^100 km s~^ for millisecond radio pulsars that have no proper 
motion measurements. 

The intrinsic P — P values of the underlying MSP population suggest that ~30% of MSPs will 
be born with apparent ages older than the age of the Galaxy. The true age distribution of MSPs 
with Tc > 10^^ yr is relatively well mixed as opposed to MSPs that lie on or just above the Hubble 
line. Hence, the sources that appear below the lO^'^yr line might not be among the oldest within 
the MSP sub-population. 

We predict younger ages for pulsars with Tc/t' >1 (Table 1). The majority of MSPs for which 
we predict younger ages, are the ones that are closer to the spin-up line. These MSPs have a 
significant fraction of their spin-down trajectories truncated because the timescales pulsars spend 
close to the spin-up line is much shorter than MSPs with smaller spin-down rates. Some of the 
younger sources with Ti < ?' < Tc are PSRs J0218-F4232, J0737-3039A, J1023+0038, B1534-M2, 
B1913-M6, and B1937-F21. 

For MSPs with no proper motion measurements. Table 2 shows the potential biases for an 
assumed 'yt=100 km s~^. In this category, PSR J1841-I-0130 may have the strongest age corruption 
with Tc/t' ~5.85. 

There are two sources in particular that have been of considerable interest: 



PSR J1012+5307: The range s of cooling ages for possibly the best studies example i s PSR 
J1012-F5307 with r^^ ~o.3-7 Gyr (jLorimer et al. Ill995l : [Alberts et all 1 19961 : burderi et a"l]ll996l l. 
We derive an age of r' ~ 6.25 Gyr (Table 1) for PSR J1012-I-5307 which is consistent with r^. 
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This implies that the true age has to be < 6.25 Gyr. One cannot exclude younger ages by 
either the Shklovskii corrected (unbiased) characteristic age (r^) or MSP age (r') approach as 
PSR J1012+5307 may have been born with a Pq very close to its currently observed period (Pq — P) 
due to low accretion rates experienced during the LMXB phase (see 

PSR B1257+12: The cumulative correction to the age of PSR B1257+12 is among the most 
significant of MSPs that have proper motion measurements {P /P' ~16.6). We infer an upper limit 
of Tt < t' = 1.42 X 10^° yr for the age. This implies that the age of PSR B1257+12 cannot be 
constrained solely from its spin-down history. 

In general, when interpreting ages inferred from spin-down histories for single MSPs, one 
has to consider the possibility that an evolutionary process which produces MSPs without a binary 
companion may affect the spin-down evolution. For instance, an encounter by which the companion 
is ejected will perturb the spin-down rate of the compact primary. In the exceptional case of PSR 
B1257+12, the process that led to the formation of the planets may have affected the spin-down 
evolution. 



4.2. Braking Index 



We tested whether alternative energy loss mechanisms other than pure dipole braking (n=3) 
are required to reconcile for the MSP age distribution. While more efficient processes (n> 3) 
may also contribute to the downward age contamination, it would seem unnecessary to invoke 
higher braking indices to account for ages that appear older than the Galaxy. We do not rule out 
that gravitational wave radiation may expedite spin-down during the very early stages after re - 



birth before magnetic stability sets in (jLindblom 



1995 



Lindblom fc Mendelllll995l : lBildstenlll998l l. 



However, for MSPs following standard spin-up, the contribution of gravitational wave radiation to 
age contamination may not to exceed the offset between the blue-to-red r age lines in Figure [TJ We 
explicitly show that lower preferred accretion rates during the active accretion phase can produce 
the paradoxically older appearing MSPs (Figures [2] and [3|) . 

Based on the P — P characteristics of MSPs, we find no compelling evidence that energy 
loss has been dominantly driven by multipole or gravitational wave radiation during a significant 
portion of the lifetime of these sources. 



5. Conclusions 

We have implemented constraints arising from the spin-up process and a limiting maximum 
spin limit into the standard method to obtain a more realistic age (r) estimate for MSPs. There 
are a range of ramifications that follow: 
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Age distribution: The unbiased characteristic ages are only upper hmits to the true age. The 
new age estimate gives tighter upper Hmits and hence is closer to the true age (rj <t'< t^). This 
flattens and shifts the age distribution toward younger ages while the age corruption scrambles the 
positions on the P—P plane quite efficiently. We predict that a significant fraction of MSPs are born 
with apparent older ages. The true age distribution of MSPs does not appear to peak at lO^'^yr 
as sharply as expected for a sub-population recycled from a first generati on of pulsars with feature s 



indicative of a population that is already old, at least in a dynamic sense (iHansen &: Phinneylll997l ). 
MSPs that appear older than the Galaxy can be reconciled with very low (m <C M^dd) progenitor 
accretion rates experienced during the latest phases of the LMXB evolution. We expect -^30% of 
the population to be born with > 10^'^ yr. 

Age corruption: There are two sources of age corruption: (1) the previously known "age bias", 
which appears to be more prominent than previously predicted and (2) "age contamination" which 
is driven by lower progenitor accretion rates. Age contamination, which effectively disguises young 
MSPs as old ones, is not correctable in the absence of additional constraints that may give us 
insight into the details of individual prior accretion histories. On the other hand, the correctable 
age bias, manifests itself as reverse contamination and will disguise old MSPs as younger sources. 
The downward contamination will remain as the main source of confusion with regards to MSP 
ages. We expect to overestimate the true age of MSPs by more than a factor of 2 for ~30% of the 
population. As a consequence, the birth and merger rates of NS-NS systems based on Tc are most 
likely underestimates which therefore will have ramifications for potential LIGO sources. 

Braking indices: The millisecond radio pulsar demographics is consistent with the canonical 
spin-down model (n=3). The challenge to disentangle possibly mixed sub-populations of MSPs that 
may have experienced dissimilar energy loss histories (n<3 or n=5) is mainly due to the paucity 
of sources. Therefore, an early short phase when MSP energy loss is dominated by gravitational 
wave or multipole emission remains as a potentially contributing source of age contamination. 



The research presented here has mad e extensive use of the 2009 August version of the ATNF 
Pulsar Catalogue (iManchester et al.ll2005l ). BK thanks Athanasios Kottas for long discussions on 
Bayesian statistics which have seeded the idea and subsequently given birth to the proper param- 
eterization of millisecond pulsar evolution. We thank the anonymous referee for useful comments. 
The authors acknowledge NASA and NSF grants AST-0506453. 
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Table 1. Ages of millisecond pulsars with proper motion measurements. 



Pulsar 


Tc [Gyr] 


< [Gyr] 


r [Gyr] 


r' [Gyr] 


Tc/r' 


J0030+0451 


7.56 


7.65 


7.24 


7.32 


1.03 


J0034-0534 


6.00 


55.71 


4.29 


39.90 


0.15 


J0218+4232 


0.48 


0.48 


O.+0.04 

'J-'^^-o.og 




1 07+0.48 
1-'J'-0.14 


J0437-4715 


1.59 


6.12 


1.47 


5.94 


0.27 


J0610-2100 


4.93 


17.92 


4.60 


16.72 


0.30 


J0613-0200 


5.06 


5.28 


4.52 


4.72 


1.07 


J0621+1002 


9.67 


10.01 


9.56 


9.91 


0.98 


J0711-6830 


5.80 


10.46 


5.61 


10.11 


0.57 


J0737-3039A 


0.20 


0.20 


'-'■^*-0.04 


Q 14+0.02 

^•1^-0.04 


, .0+0.57 
i-^''-a.i8 


J0751+1807 


7.08 


7.25 


6.49 


6.66 


1.06 


J1012+5307 


4.87 


6.48 


4.69 


6.25 


0.78 


J1023+0038 


2.23 


2.50 


1.45 


1.62 


1.37 


J1045-4509 


6.77 


10.91 


6.63 


10.71 


0.63 


B1257+12 


0.86 


14.58 


0.75 


14.21 


0.06 


J1453+1902 


7.91 


8.46 


7.68 


8.21 


0.96 


J1455-3330 


5.21 


8.07 


5.07 


7.93 


0.66 


J1518+4904 


23.84 


29.34 


23.74 


29.23 


0.82 


B1534+12 


0.25 


0.25 


n 1 q+001 


on+O Ol 
^■^^-0.04 


1 94+0-32 
1-^^-0.05 


J1600-3053 


6.00 


6.76 


5.54 


6.24 


0.96 


J1603-7202 


14.98 


17.94 


14.85 


17.81 


0.84 


J1640+2224 


17.71 


30.59 


15.94 


27.53 


0.64 


J1643-1224 


3.96 


5.04 


3.77 


4.81 


0.82 


J1709+2313 


20.21 


49.45 


19.27 


47.15 


0.43 


J1713+0747 


8.49 


9.01 


8.08 


8.58 


0.99 


J1738+0333 


3.85 


4.07 


3.71 


3.93 


0.98 


J1744-1134 


7.24 


9.36 


6.80 


8.80 


0.82 


B1855+09 


4.80 


4.92 


4.63 


4.75 


1.01 


J1909-3744 


3.34 


17.08 


2.95 


15.12 


0.22 


T1 01 1 111/1 


4.U0 


y. io 


0. / 4 


0.44 


(\ AQ 


B1913+16 


0.11 


0.11 


07+0 01 
'-'■^'-0.03 


07+0-01 
^■l"-0.03 


1 66+1-0" 
i.DD_Q 29 


B1937+21 


0.24 


0.24 


in+0 04 


in+o-04 

^•1^-0. 09 


r, 07+21.63 
^•■^'-0.66 


J1944+0907 


4.80 


8.59 


4.63 


8.27 


0.58 


B1953+29 


3.27 


3.41 


3.14 


3.27 


1.00 


B1957+20 


1.51 


2.23 


0.92 


1.37 


1.10 


J2019+2425 


8.88 


24.34 


8.31 


22.77 


0.39 


J2051-0827 


5.63 


5.81 


5.35 


5.52 


1.02 


J2124-3358 


3.79 


6.25 


3.64 


5.99 


0.63 


J2129-5721 


1.98 


2.09 


1.84 


1.94 


1.02 


J2145-0750 


8.53 


9.81 


8.42 


9.69 


0.88 


J2235+1506 


5.99 


9.13 


5.92 


9.04 


0.66 


J2317+1439 


22.55 


36.18 


20.65 


33.13 


0.68 


J2322+2057 


7.85 


18.49 


7.51 


17.69 


0.44 



Note. — PSRs J1024-0719, J1801-1417 and J2229+2643 have P' < for 
Dt=146, 208 and 115 km s~i respectively. Tc and r are the biased ages while 
and r' are unbiased for the effects of secular acceleration. B1937+21 is 
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Fig. 1. — Observed distribution of Galactic MSPs: black dots (•) indicate unbiased spin-down rates 
{P') for MSPs that have proper motion measurements. Vertical solid lines are the corresponding 
Shklovskii correction. MSPs with no proper motion measurements (o) are corrected for an assumed 
vt = 100 km s~^ in order to see the potential level of bias. The magenta squares are MSPs (■: PSRs 
J1022+1001, J1216-6410, J1829+2456 and J1933-6211) for which the spin-down rates are corrected 
only for an assumed vt = 50km s~^ while these MSPs will otherwise appear to be spinning-up, 
i.e. < 0. The green triangles are MSPs (a: PSRs J1024-0719, J1801-1417 and J2229-F2643) 

with proper motion measurements for which P' < and therefore are left uncorrected. The red 
and blue lines are the MSP age (r) lines for braking indices n=3 (- - ' ) and 5 (, • ' ). The diagonal 
solid line is the spin- up line P ~ ttiPq''^ for ?h = MecLcL- The diagonal dotted lines ("•• ) are the 
inferred dipole field strengths. {~--~ ) is the critical magnetic field below which MSPs can be 
spun up to the mass shedding limit Pgh —1 ms. 
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Fig. 2. — Expected true age distribution of the underlying MSP population: color represents the 
true age (tj). Downward arrows {],) are the neutron star spin frequencies measured in LMXBs 
which are used as progenitor seeds to reconstruct the synthetic population. 
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Fig. 3. — P — P and underlying age distribution of MSPs (a) before and (b) after correcting for 
effects of secular acceleration: color represents the true age {jt) and the diagonal dashed lines 
are the inferred characteristic age (tc) lines. The progenitor seeds used to produce the synthetic 
population are randomly sampled from the period distribution of observed millisecond X-ray pulsars 
(downward arrows), (a) The transverse velocities and distances are chosen to be consistent with 
the observed millisecond radio pulsar population in order to see the potential bias in the P — P 
demographics, (b) The P — P distribution for the same sample MSP population after correcting 
for effects of secular acceleration. The process that disguises older MSPs as younger sources can be 
reversed by properly correcting (unbiasing) for the Shklovskii effect, whereas the population will 
still harbor young MSPs with apparent older ages. This (downward) age contamination is driven by 
the lower mass accretion rates experienced during the LMXB phase. We predict that a significant 
fraction of millisecond radio pulsars (~30%) will be born with characteristic ages older than the 
age of the Galaxy. 
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Fig. 4. — Effects of secular acceleration on the MSP age distribution (see also Figures [3]): ages 
inferred from measured P and P values will overestimate or underestimate the true age (dashed 
line). Unbiased characteristic ages (r^) reflect mere overestimates (solid line) where -^30% of MSPs 
will have a overestimating the true age by more than a factor of 2. 
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likely to be significantly younger than its characteristic age. Associated errors 
that are more than 5% of the most likely value are tabulated. Calculations are 
performed in double precision before rounding. 
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Table 2. Ages of millisecond pulsars with no proper motion measurements. 



Pulsar 


Tc [Gyr] 


[Gyr] 


r [Gyr] 


r' [Gyr] 
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l.Uz 


Ti A on c;ftoc; 
J 14zU — ODZo 


o.Ui 


1 1 fiT 


v ni 
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4.46 
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4.36 
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44 


u.^o 


Tin 
1 . ±u 


TI 9940 
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i.yo 
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J1756-2251 


0.44 


0.45 


0.38 


0.38 


1.16 


J1757-5322 


5.34 


7.30 


5.21 


7.16 


0.75 


J1802-2124 


2.78 


2.95 


2.68 


2.84 


0.98 


J1804-2717 


3.62 


4.59 


3.50 


4.46 


0.81 


J1810-2005 


3.44 


3.66 


3.36 


3.57 


0.96 


J1841+0130 


0.058 
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3.41 
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2.41 


1.27 
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7.33 


10.65 


6.89 


10.01 


0.73 
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1.81 


1.85 


1.42 


1.45 


1.25 


J1904+0412 


10.24 


12.40 


10.16 


12.32 


0.83 


J 1905+0400 


12.34 


33.12 


11.47 


30.81 


0.40 


J1910+1256 


8.08 


11.27 


7.76 


10.81 


0.75 


J1911+1347 


4.29 


5.24 


4.09 


4.99 


0.86 
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5.05 


6.69 


4.91 


6.55 


0.77 


J2010-1323 


17.17 


185.03 


16.54 


178.24 


0.10 


J2033+17 


8.57 


14.86 


8.33 


14.44 


0.59 



Note. — We assume Dt = 100km s ^ for the Shklovskii correction in order 
to see the potential bias. PSRs J1022+1001, J1216-6410, J1829+2456 and 
J1933— 6211 have P' < for i)t=100km s^^. Tc and t are the biased ages while 

and t' are unbiased for the effects of secular acceleration. PSR J1841+0130 
is likely to be significantly younger than its characteristic age. Associated errors 
that are more than 5% of the most likely value are tabulated. Calculations are 
performed in double precision before rounding. 



